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I.   INTRODUCTION 

A.   GENERAL 

The  helicopter  design  process,  the  subject  of  numerous 
articles  and  studies  is  an  evolving  discipline  that  borders 
on  being  an  art.   A  successful  design  must  balance  the 
user's  needs  and  desires  against  practical  capabilities. 

With  the  introduction  of  composite  materials  and  new 
technologies,  principally  in  rotor  and  engine  performance, 
significant  advances  have  been  made  in  helicopter  capabil- 
ities.  In  some  instances,  the  performances  of  hybrid 
helicopter  designs  rivals  that  of  a  similarly  sized  conven- 
tional aircraft.   For  example,  the  YVX,  a  joint  Boeing-Bell 
venture,  will  have  the  hover  and  low  speed  capabilities 
of  a  helicopter  while  being  able  to  cruise  at  300  knots. 

Viable  commercial  and  military  helicopter  designs  are 
only  thirty  years  old.   The  first  major  use  of  helicopters 
occurred  during  the  Korean  conflict.   To  put  this  in 
perspective,  the  first  large  scale  use  of  conventional  type 
aircraft  was  in  World  War  I. 

Helicopter  design  can  proceed  on  a  number  of  different 
levels,  ranging  from  comprehensive  computer  design  programs 
to  preliminary  analysis  using  simplifications  of  the  basic 
performance  equations.   Each  has  its  merit  and  place. 
Computer-aided  design  provides  a  great  deal  of  data. 
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Generally,  these  programs  integrate  aircraft  configuration 
sizing,  performance  and  weight  calculations  in  an  iterative 
process.   An  example  of  a  computer  design  program  for  heli- 
copters is  the  Helicopter  Sizing  and  Performance  Computer 
Program  [HESCOMP] ,  orginally  developed  by  Boeing- Vertol 
for  NASA.   This  program  is  currently  used  as  a  wide  number 
of  institutions  conducting  studies  in  helicopter  design. 

On  the  opposite  end  of  the  spectrum  would  be  sensitivity 
design  studies  using  the  performance  equations.   Surprisingly 
accurate  simplications  of  these  equations  can  be  made.   This 
provides  the  designer  with  an  excellent  method  for  doing 
first  cut  preliminary  helicopter  sizing  at  a  low  cost. 

B.   OBJECTIVE 

This  report  is  an  investigation  of  several  of  the 
methods  employed  in  the  preliminary  design  of  a  helicopter. 
Conceptually,  the  report  can  be  divided  into  three  parts. 
In  the  first  section,  a  sensitivity  analysis  of  the  basic 
performance  equations  is  performed.   The  purpose  here  is  to 
ascertain  where  reasonable  simplications  can  be  made 
that  do  not  seriously  degrade  the  accuracy  of  the  result. 

In  the  second  section  a  graphical  method  of  doing 
parametric  design  studies,  known  as  Carpet  Plots,  is 
developed.   This  method  allows  the  user  to  formulate  a 
graphical  solution  matrix  to  meet  the  design  criteria 
specified  for  the  helicopter.   Carpet  Plots  are 
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particularly  instructive  since  they  give  visual  insight 
into  the  interplay  of  the  various  design  parameters. 

In  the  last  section,  an  overview  of  HESCOMP  is  given 
Programs  are  developed  which  enable  a  person  to  access 
HESCOMP  on  the  Naval  Postgraduate  School  Main  Frame  IBM 
system. 
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II.   SENSITIVITY  ANALYSES  OF  BASIC  HELICOPTER  EQUATIONS 

A.  DESCRIPTION  OF  PROBLEM 

In  preliminary  helicopter  design,  there  are  a  number  of 
instances  where  a  quick  first  cut  analysis  would  be 
extremely  helpful.   This  is  especially  true  in  determining 
the  preliminary  size  of  the  helicopter  required  to  meet 
the  specifications. 

Historically,  there  are  a  number  of  variables  in  the 
performance  equations  of  helicopters  which  may  be  treated 
as  constants.   This  may  allow  for  significant  simplifica- 
tions and  aid  in  the  preliminary  design  process. 

In  this  section,  a  sensitivity  analysis  of  the  per- 
formance equations  is  done.   In  a  sensitivity  analysis, 
each  parameter  [or  variable]  is  varied  in  order  to  deter- 
mine its  effect  on  the  equation.   Variables  which  are 
shown  to  have  little  effect  may  be  treated  as  constants  and 
the  equation  simplified  accordingly. 

B.  SOLIDITY 

Solidity,   a  ,  is  the  fraction  of  the  disk  area  that  is 
composed  of  blades.   It  is  a  function  of   b  ,  the  number 
of  blades,  of  a  constant  cord,  c  ,  at  a  radius,  R: 


a  =  5|  (2.1) 

TT  R 
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C.   DISK  LOADING: 

Disk  loading  is  defined  as  the  ratio  of  the  weight  to 
the  total  area  of  the  rotor  disk. 


m    WEIGHT 
UL    ~   AREA 

(2.2) 

-  £ j    [lo/ft  ] 

ttR 


D.   POWER  LOADING 

Power  loading  is  the  ratio  of  weight  to  input  power. 

PL  =  p^  [lb/hp]  (2.3) 

in 

In  a  hover,  thrust  equals  weight;  this  allows  us  to 
rewrite  the  power  loading  for  the  hover  condition  as 


PT    _L-     ROTOR  THRUST    nh/h    i        r?  d1 
P.     ROTOR  HORSEPOWER  L1D/nP-l        U-^J 
in 


E.   COEFFICIENT  OF  THRUST  AND  POWER 

The  coefficient  of  thrust,   L  ,  is  a  non-dimensional 
coefficient  which  facilitates  computations  and  comparisons 


C   =  — ^  =  —2-1 2  (2.5) 

ApV-jT    ttR  p(ftR) 


Similarly,  a  coefficient  of  power,   C   ,  has  been 
established  as : 
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— ; z         (2-6) 

ApVT"    ttR"p(2R) 


No  significant  simplifications  can  be  made  to  either  of 
these  coefficients.   However,  it  should  be  observed  that 
the  coefficient  of  thrust  is  inversely  proportional  to  the 
square  of  the  rotor  tip  velocity,  while  the  coefficient  of 
power  is  inversely  proportional  to  the  cube. 

Assuming  all  other  factors  being  equal,  increasing  the 
rotor  tip  velocity  from  600  fps  to  700  fps  [an  increase  of 
16.7  percent]  will  have  the  following  result  on  these 
coefficients . 


C 


T        2 
ApV^T 


(2.5) 


Ap(1.167)Z 


T 


Ap  (1.361) 


The  coefficient  of  thrust  is  reduced  by  26.9  percent 
Similarly,  for  the  coefficient  of  power: 
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c.  -     -  ~ 

P    .  „  3 

ApVT 


Ap(1.167) 


3  (2.6) 


Ap(1.589) 
The  coefficient  of  power  is  reduced  by  37.1  percent. 

F.   HOVER  POWER 

The  total  power  in  a  hover  is  made  up  of  two  terms, 
profile  power,   P   ,  and  induced  power,  P.  . 

Utilizing  black  element  theory  the  profile tpower 
required  to  hover  can  be  expressed  as : 

The  induced  power  predicted  by  momentum  theory  is: 


P.  =  V.   T 
i     in 


T5/2 


(2.8) 


/-y       d2 

ZTTpR 


The  total  power  required  to  hover  is 


PT  -  ?i  ♦  P0  (2.9) 
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T3/2     1     -  3 

PT  =      ^  +  ±   ar  CdQ  p  A(flR)*         (2.10) 

v  ZnpR2 


Donald  M.  Layton  in  Helicopter  Performance,  [Ref.  1], 
found  that  for  the  optimum  hover  power,  the  induced  power 
is  equal  to  twice  the  profile  power.   The  analysis  was 
performed  in  the  following  manner. 

By  assuming  constant  weight,  density,  solidity,  and  an 
average  profile  drag  coefficient,  as  well  as  a  fixed 
rotational  velocity,  equation  (2.10)  reduces  to 


P  =  T~   +  C2r2  (2'11} 


where   C-,   and   C?   are  constants. 

As  equation  (2.12)  shows,  profile  power  increases  as 
the  square  of  the  blade  radius  while  the  induced  power 
decreases  with  increasing  blade  radius. 

The  optimum  hover  power  with  respect  to  rotor  radius 
can  be  determined  by  taking  the  differential  and  setting  it 
equal  to  zero. 

m-    °  "  "  |+  2  C2R  <2-12A) 

or  —  =  2  C2R2  (2.12B) 

which  implies  Pi  =  2  Po  (2.12C) 
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G.   HELICOPTER  SIZING 

A  simplified  relationship  between  the  total  power 
required,  gross  weight  and  rotor  radius  can  be  developed  in 
the  following  manner. 

The  total  power  required  to  hover  equation  for  the  main 
rotor  was  developed  in  the  preceding  section  and  is 
repeated  here  for  clarity. 

PT  -  Pi  +  P0  (2.9) 

5/2 
PT  ■  —  •  5  +  5»r  Udo  p*  vtip3  r2       <2-10) 

/2tt  p  ^ 

In  a  hover,  thrust  equals  weight.   Solving  equation 
(2.11)  for  weight  one  obtains: 


W3/2  =  [PT  -  -g-  a  Cdo  pit  VT^  R2]  /pA       (2.15) 


This  equation  may  be  further  simplified  if  it  is 
assumed  that  the  density,  average  profile  drag  coefficient 
and  tip  velocity  are  constants;  these  are  reasonable  assump 
tions.   Historically,  the  average  profile  drag  coefficient 
of  a  helicopter  has  been  approximately  0.01.   The  operating 
environment  of  today's  helicopters,  especially  military, 
is  below  5,000  feet  agl.   This  allows  for  the  use  of  the 
standard  sea  level  value  for  density  with  little  error. 
Primarily,  due  to  tip  mach  effects,  the  upper  limit  on  the 
rotor  tip  velocity  is  in  the  range  of  700  f ps . 
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The  resulting  equation  with  these  assumptions  incor- 
porated into  a  constant,  K  ,  is: 

W  =  [47.527  PT  R  -  JH      bc]2/3         (2.13) 

Equation  (2.13)  can  be  further  reduced  when  the  order 
of  magnitude  of  the  two  terms  is  considered. 

47.527  PT  R  >>  K,  be 

Thus  , 

W  %  [47.527  PT  R]2/3  (2.14) 

To  determine  how  accurate  this  simplification  is,  the 
equation  is  used  to  approximate  the  total  weight  of  a 
number  of  helicopters  for  which  the  parameters  are  available, 
As  Table  2.1  indicates,  the  weight  approximation  formula 
yields  values  within  six  percent  of  the  actual  total  weight 
of  these  helicopters. 

H.   FIGURE  OF  MERIT 

A  figure  of  merit,  FM  ,  has  been  defined  for  the 
helicopter  as  the  ratio  of  the  ideal  rotor  induced  power  to 
the  actual  power  required  to  hover,  with  non-uniform  induced 
velocity,  tip  losses  and  profile  drag  power. 
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TABLE  2.1 


HELICOPTER  WEIGHT  COMPARISON 


HELICOPTER 

TOTAL  GROSS 

CALCULATED 

PERCENT  OF 

WEIGHT 

GROSS  WEIGHT 

ACTUAL 

(1000  lbs) 

(1000  lbs) 

GROSS  WEIGHT 

AH-64 

14.66 

14.69 

101% 

UH-1N 

14.20 

13.74 

97% 

H-3H 

21.00 

20.63 

98% 

S76 

10.00 

9.90 

99% 

UH-6DA 

20.25 

19.53 

95% 

• 

H-543 

42.00 

42.00 

100% 

H-53D 

42.00 

41.00 

98% 

H-53E 

73.50 

69.00 

84% 
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In  a  hover,  the  figure  of  merit  may  be  written  as: 


FM  = 


rr 


PL 


DL 


/p" 


(2.15) 


CT- 


/2C 


The  figure  of  merit  is  customarilty  plotted  against  the 
quantity   CT/a  .   According  to  Zalesch  [Ref .  2] ,   CT/a  ,  is 
proportional  to  the  average  blade  angle  of  attack  and  can 
be  used  as  a  measure  of  rotor  efficiency.   The  curve  in 
Figure  2.1  is  based  on  data  from  Reference  2  for  a  typical 
tail  rotor  helicopter. 

Main  Rotor  Hover  Performance 


0.00 


0.04 


008 


0.12 


Blade  Loading 


0.18 


Figure  2.1.   FM   Versus  Blade  Loading   CT/a 
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Previous  studies  have  shown  that  a  figure  of  merit 
between  0.70  and  0.80  is  considered  average . [Ref .  5] 
If  the  induced  power  is  between  70  and  80  percent  of  the 
total  power,  the  figure  of  merit  will  be  approximately  0.75. 

With  the  figure  of  merit  limited  to  values  between 
0.70  and  0.80,  the  following  simplification  can  be  made, 
assuming  the  hover  condition  of  thrust  equaling  weight  and 
standard  sea  level  conditions: 


FM  = 


IV 


3/2 


67.214  PTR 


(2.16) 


For  Navy  helicopter  design,  the  rotor  radius  has  been 
limited  by  flight  deck  spotting  constraints  to  less  than 
30  feet;  the  exception  to  this  is  the  H-3,  R  =  51  feet  and 
the  H-53,  R  =  56  to  58  feet  [depending  on  the  model]. 
However,  these  two  helicopters  work  almost  exclusively  from 
large  air  dedicated  ships  such  as  the  LPH,  LHA  and  CV. 

If  the  small  deck  operating  assumption  is  made, 
equation  (2.16)  can  be  further  simplified  to  [assuming 
R  =  28  feet] : 


P  = 


W 


5/2 


1881.98  FM 


(2.17) 


An  FM  of  0.80  will  yield  a   P   to   W  relationship  of 


W 


5/2 


T    1505.58 


(2.18) 
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while  an   M   of  0.70  yields  a  relationship 


w3/2 

P  "  1517.39  f-2  '  19^ 


If  equation  (2.17)  is  solved  utilizing  the  approximate 
weight  relationship  developed  earlier  of 


W3/2  =  47.527  PTR  (2.  14) 


a  value  for  the  figure  of  merit  of  0.707  is  obtained.   This 
is  within  the  historical  range  of  values. 

I.   TAIL  ROTOR  SIZING 

A  historical  analysis  of  typical  helicopters  [Ref.  3), 
shows  the  following  empirical  relationship  for  the  tail 
rotor  radius 


PT  *  1.3  [x^]1/2  [ft]  (2.20) 


when  comparing  the  results  of  this  equation  with  actual 
tail  rotor  radius  data,  it  was  found  that  if  a  multipli- 
cation factor  of  1.2  is  used  vice  1.5  a  better  approximation 
is  obtained.   The  results  are  tabulated  in  Table  2.2. 

J.   FORWARD  FLIGHT  POWER  CONSIDERATIONS 

The  total  power  in  forward  flight  consists  of  induced, 
profile  and  parasite  power.  If  the  helicopter  is  a  single 
rotor  vehicle,  the  tail  rotor  power  should  be  taken  into 


23 


HELICOPTER 

AH- 

■64 

UH- 

-IN 

SH- 

■3H 

s-: 

76 

UH- 

■6  0A 

CH- 

-53D 

CH- 

-53E 

TABLE  2.2 


TAIL  ROTOR  SIZING 


ACTUAL  TAIL  ROTOR      APPROXIMATION  [FT] 
RADIUS  [FT]  [2.20]       [2.21] 


4.6  4.98  4.59 

4.3  4.90  4.52 

5.3  5.95  5.5 

4.0  4.11  3.79 

5.5  5.85  5.4 

8.0  8.42  7.78 

10.0  11.15  10.29 
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account,  as  well  as  all  mechanical  losses  [transmission, 
etc.]  for  accurate  calculations.   However,  a  reasonable 
approximation  can  be  obtained  by  considering  only  the  main 
rotor  and  increasing  this  power  figure  by  several  percent 
to  account  for  these  losses. 

Figure  2.2  is  a  plot  of  the  induced,  profile,  parasite 
and  total  power  curves  for  typical  tail  rotor  helicopter. 

MAIN  ROTOR  POWER 
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Figure  2.2.   Power  Required  Versus  Forward  Velocity 

The  induced  power  drops  off  rapidly  with  increasing 
forward- velocity ,  whereas  the  parasite  power  increases 
rapidly. 

Parasite  power  is  the  power  required  to  overcome  the 
drag  forces  created  by  the  aircraft's  geometry.   These  drag 
forces  are  due  to  pressure  drag  and  skin  friction. 
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Parasite  drag  is  extremely  sensitive  to  the  helicopter's 
loading.   It  is  generally  a  minimum  for  forward  flight  and 
increases  for  sideways  flight.   Helicopters  are  generally 
streamlined  for  forward  flight  and  the  flat  plate  area  is 
a  minimum  in  this  direction.   The  equation  for  the  parasite 
power  is: 

Pp  =  \   P  V£3  ff  (2.2  1) 

The  parasite  power  is  a  function  of  the  cube  of  the 
forward  velocity.   As  such,  with  the  advent  of  high  speed 
helicopters  a  great  deal  of  consideration  has  been  placed 
on  streamlining  the  geometric  shape  in  order  to  reduce  this 
power  requirement. 

Blade  element  theory  is  commonly  used  to  develop  the 
profile  power  equation  for  forward  flight.  An  excellent 
development  of  this  equation  is  given  in  Reference  1. 

The  profile  power  equation  in  forward  flight  is: 

Pof  =  i  a  CdQ  p  A  VT3  [1  +  4.5  y2]         (2.22) 

Equation  (2.23)  is  a  function  primarily  of  the  main 
rotor  geometry.   The  variable  with  the  most  significance 
is  the  rotor  tip  velocity;  increasing  the  tip  velocity  from 
600  to  700  fps  results  in  a  58.8  percent  increase  in 
profile  power  [assuming  other  factors  are  constant] . 
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The  induced  power  is  a  function  of  the  induced  velocity. 
In  a  hover,  the  total  flow  through  the  rotor  system  is 
induced.   As  the  forward  velocity  increases,  the  mass  flow 
rate  through  the  rotor  disc  increases  due  to  the  forward 
translation  of  the  helicopter.   This  reduces  the  induced 
velocity. 

The  equation  for  the  induced  power  requirements  at  all 
forward  velocities  is  : 


P  =  T  .  V.* 
it 


(2.2  3) 


where 


V 


it 


vf2/v2 


?    2 
V£V2VZ 


1/2 


T } 

+  u      .v 


(2.25a) 


At  high  forward  velocities,  the  induced  power  required 
can  be  approximated  as: 


W 

P.  =  wv.+  =  "v 

i      it    2pAVr 


(2.24) 


The  total  power  for  forward  flight  is  the  sum  of  the 
induced,  profile  and  parasite  powers. 


PT  "  Pi  +  Po  +  ?p 


(2.25) 
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PT    =    T.Vit    +    ^  a    CdQ    p    A   VT3[1    +    4.3   y2] 


(2.25a) 


*  \  P    ££  vf3 


At  high  forward  velocities,  equation  (2.23)  can  be 
substituted  into  equation  (2.25),  resulting  in: 


W 


T    2pAVr    8  u  "do 


♦  \  P  ff  vf3 


1  3         Vf 

i  a  C,   p  A  V/[1  +  4.3   ±] 


T 


P,R- 


I£  one  makes  the  following  assumptions 


W  =  const    C,   =  const 


(2.26) 


p  =  const 
VT  =  const 


a  =  const 


Equation  (2.26)  reduces  to 


Kl     2   2 
PT  =  -4  +  YT    R   +  P 

T    R2  p 


(2.27) 


The  derivative  of  equation  (2.27)  with  respect  to  radius  is 


dP 


T 


2K- 


dR 


R 


t   +    2    K2  R 


(2.28) 
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Setting  this  equal  to  zero,  one  obtains 


2K 

--±  +  2  K?  R  =  0  (2.28a) 

R        Z 


R       2K1 

7  *  [-  -V  +  2  K?  R]  =  0  (2.28b) 

R-5        Z 


Kl        2 

-4  =  K?  IT  (2.28c) 

R^     Z 


P.  =  P  (2.28d) 

10 


This  defines  point  of  minimum  total  power  required  for 
VMAX  range.   This  corroborates  with  the  res.ults  obtained  by 
Waldo  Carmona  [Ref.  4]. 

If  the  total  power  required  is  differentiated  with 
respect  to  forward  velocity  and  is  set  equal  to  zero,  it  can 
be  seen  that 


P.  =  5  P  (2.29) 

10 


or 


w2      3p£Vf 

_^ =  L_  (2    30) 
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Solving  this  equation  for  velocity  results  in 


vf  ■ 


[( 


W     A 


1/21  1/2 


A   3f 


ft/sec 


(2.31) 


According  to  Carmona  [Ref .  4] ,  this  corresponds  to  the 
best  endurance  velocity. 

K.   DENSITY  EFFECTS  ON  TOTAL  POWER 

The  effect  of  density  on  the  total  power  required  in 
forward  flight  is  as  follows: 

The  general  operating  altitudes  of  a  helicopter  are 
below  10,000  feet.   The  corresponding  ICAO  STANDARD 
ATMOSPHERE  range  for  density  is 

p  =  0.0025769  [lb  sec2/ft4]   SSL 

p  =  0.0017553  [lb  sec2/ft4]   at  10,000  feet 

p/pSSL   varies  from  1  to  .7385. 

The  effect  on  the  components  of   PT   are  as  follows: 
Induced  Power: 

1/p/pSSL  =>  l  to 


7385 


This  translates  to  a  35  percent  increase  in  the  induced 
power . 
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Parasite  and  Profile  Power: 

Both  parasite  and  profile  powers  are  directly  propor- 
tional to  the  density  ratio.   Therefore,  as  you  go   up  in 

altitude  both   P    and   P    are  reduced. 

o        p 
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III.   CARPET  PLOT  DESIGN  STUDY 

A.   DESCRIPTION  OF  PROBLEM 

Preliminary  helicopter  design  involves  one  with  a  wide 
range  of  choices.   For  any  given  payload  and  performance 
specifications,  there  a  number  of  helicopter  designs  that 
satisfy  the  requirements.   The  problem  in  the  preliminary 
design  process  is  narrowing  these  possibilities  and 
selecting  the  design  which  will  provide  the  best  helicopter 
for  the  mission. 

Obviously,  the  operating  environmental  constraints  help 
to  define  the  basic  configuration.   These  constraints  are 
usually  specified  in  the  Request  for  Proposal  [RFP],  in  the 
case  of  a  military  helicopter.   For  example,  typical 
constraints  placed  on  the  design  of  a  Navy  helicopter  are 
the  size  of  the  ship  deck  and  hangar  from  which  it  will  be 
operating,  the  requirement  for  a  blade  fold  system,  dual 
engine  configuration  and  IFR  capability. 

Even  with  these  design  constraints,  there  is  still 
a  great  deal  of  leeway.   In  order  to  insure  that  the  best 
helicopter  design  is  selected,  an  appropriate  number  of 
solutions  satisfying  the  specifications  should  be  inves- 
tigated.  Since  each  solution  is  generally  characterized 
by  a  different  combination  of  design  parameters,  the 
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selection,  according  to  Greenfield  [Ref.  5],  can  best  be 
made  through  a  parametric  study  which  allows  for  the 
optimization  of  many  design  parameters. 

One  method  of  parametric  analysis  used  is  Carpet  Plots. 
This  method  is  based  on  the  simultaneous  graphical  solution 
of  the  weight  and  hover  performance  equations.   To  this 
solution  set  is  added  to  the  environmental  constraints  to 
the  helicopters  size.   This  effectively  brackets  the  area 
of  acceptable  design  solutions. 

This  method  assumes  that  minimum  gross  weight  is  the 
criterion  by  which  the  best  [or  optimum]  design  parameters 
are  selected. 

B.   ASSUMPTIONS 

1.  Airfoil  used  is  a  derivative  of  the  NACA  0012 

with  the  following  mean  approximate  values  from  Reference  5 
a  =  slope  of  airfoil  section  lift  curve,   dC  /da  , 

per  rad. 
a  =5.73 

6   =  blade  section  drag  coefficient 
5Q  =  .009 

62  =  .3 

2.  a)   The  tail  rotor  radius  is  assumed  to  be  .16 
times  the  main  rotor  radius  [Ref.  5]. 
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b)   The  distance  between  the  rotors,  or  tail  rotor 
moment  arm,   £TR   is  1.19R  [Ref.  5]  .   These  ratios  reflect 
the  values  of  maximum  rotor  diameter  and  overall  length 
specified  as  size  limitations. 

3.   B  =  .97.   Historical  approximation  [Ref.  7]. 

C.   METHODOLOGY 

In  order  to  properly  develop  the  weight  and  performance 
equations  required  for  a  carpet  plot  design  study,  the 
payload  and  performance  specifications  of  the  helicopter 
are  needed.   This  data  is  used  to  tailor  the  equations  for 
the  design. 

The  equations  will  be  developed  here  for  a  four-place 
light  helicopter.   The  equation  development  procedure  is 
applicable  to  other  size  helicopters;  the  development  for  a 
medium  helicopter,  20,000  lb  weight  class,  is  to  be  found 
in  Appendix  B. 

The  following  specification  requirements  which  are 
similar  to  those  in  Reference  5  will  apply  to  this  design: 

1.  The  rotor  diameter  should  be  less  than  35.2  feet. 

2.  The  overall  length  should  be  less  than  41.4  feet. 

3.  The  gross  weight  of  the  helicopter  should  not 
exceed  2,450  lbs. 

4.  The  helicopter  should  be  capable  of  hovering,  out 
of  ground  effect  at  6,000  feet  with  an  ambient  air  tem- 
perature of  95  F. 
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5.  The  useful  load  at  hover  shall  consist  of,  as  a 
minimum,  200  lbs  for  the  pilot,  400  lbs  of  payload  and 
sufficient  fuel  to  give  the  helicopter  up  to  three  hours 
endurance  at  sea  level  conditions. 

6.  Maximum  speed  of  at  least  110  knots  using  Normal 
Rated  Power,  at  sea  level. 

7.  Total  Power  Required  at  6,000  feet  and  95  F  shall 
be  not  more  than  206. 

D.   HOVER  EQUATIONS 

1.   The  main  rotor  power  required  to  hover  out  of 
ground  effect  is 

Total  Main  Rotor  Power  [Hover]  =  Rotor  Profile  Power  +  Rotor 
Induced  Power 


1 . 1  3W 


DL 


550B/2pQ  ,/p/po 


6WY"T  p/p 


4400  C 


LPo 


C 


LRo 


0      I ap/pQ  J 


(3.1) 


At  an  altitude  of  6,000  feet  and  a  temperature  of  95   , 

p/p   =  .749395  .   Therefore,  equation  (1)  can  be  simplified 

to  : 


PT6000/95°F  "  -035479W[DL] 


1/2 


(3.2) 


+    -^1971     [10]"5(1    +    1. 
LLRo 


8  0  779    CT  n    ")W   VT 
LRo  I 
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The  tail  rotor  thrust  required  to  counterbalance  the 
main  rotor  torque  is : 


550  PTR    550  PT 


TR 


£TR  VT 


1.19V. 


(3.3) 


where   &TR  has  been  defined  as  1.19R  .   With   R„R   defined 
as  . 16R  ,  the  tail  rotor  disk  loading  can  be  written, 
using  equation  (3)  as: 


DL 


TR 


550    P 


T 


TR        ATD         1.19VT       \    77ZZ2 
iR  T   tt  (  .  16R) 


(3.4) 


550    P 


T  DL 


1.19(.0256)VT   W 


Greenfield  [Ref.  5],  in  his  development,  assumes  that 
the  tail  rotor  tip  speed  is  equal  to  the  main  rotor  tip 
speed  and  that   oTR  =  .02   and   3TR  =  .90  .   With  these 
assumptions  the  equation  for  the  tail  rotor  power  required 
to  hover  can  be  written  as : 

13/2 


=  2055.7 


TR 


Hover 


DL 

W  p/p 


1/2 


.012605  P 
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LRTR 
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T 


(3.5) 
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The  equation  for  the  tail  rotor  mean  blade  lift  coeffi 
cient  can  be  written  as 


T 


LRTR    562.5(p/p  ) 


(3.6) 


if  it  is  assumed  that  the  tail  rotor  is  designed  to  counter 
balance  a  sea  level  main  rotor  torque  equivalent  to  90 
percent  of  the  installed  power. 

Substituting  equation  (3.6)  into  equation  (3.5)  one 
obtains  the  following  expression  for  hover  tail  rotor  power 


=  2574.7 


TR6000/950 


DL 
W 


1/2 


T 


II 


V 


T 


3/2 


+  5.3134       (5.7) 


It  is  assumed  that  the  gear  losses  amount  to  3  percent 
and  that  there  is  a  1  percent  cooling  power  loss,  the  total 
brake  horsepower  required  to  hover  becomes: 


T 


P    +  P 
rTm    TTR 

96 


(3.8) 


Empirical  studies  have  shown  that  the  tail  rotor  power 
required  to  hover  can  be  approximated  by 


AC 


8  [total  horsepower  to  hover] 


This  allows  one  to  write  the  main  rotor  power  required  to 
hover  as  : 


PTm  "  C-88^PTnP 


(3.9) 
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Following  Greenfield's  [Ref.  5]  development  further, 
if  equations  (5.2)  and  (3.7)  are  substituted  in  equation 
(5.8) ,  one  obtains 


PT  =  .036757  W  /DL 

^6000/95° 


(10)"5  [1  +  1.80779  CTD^2]  W  VT   (3.10) 


'LRo 


+  2475.6 


LRo 


'DL 
'W 
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3/2 


+  5.534 


Utilizing  the  approximation  for  tail  rotor  power, 
equation  (3.9),  equation  (3.10)  can  be  solved  for  W 
(gross  weight)  as  a  function  of  variables   VT  (tip  speed) , 
DL   (rotor  disk  loading),   CfR   (rotor  mean  lift  coefficient) 


an 


d   PT   (total  power  to  hover) 


W    = 


411.51 


DL 


i/4 


V, 
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7  "1 


TJl 


1    +    K. 
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/DL 


-    K 


VT    +    K4    /DL 


(5.11) 


where 


K,    =   P 


(10) 
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t r  o  5  o  4  b  0 

K3  =  — K—  LZ-Hc) 


K      _  5695.7 

K4 k (3. lid) 


K5  "  4^  CI  *  L80779  CLRo2)  (3. lie) 

LRo 


Equation  (3.11)  has  been  programmed  in  Appendix  B 
and  solved  for  tip  speeds  from  600  to  700  cps  and   CT D 
of  .3  to  . 7 . 

Equation  (3.11)  is  one  of  the  two  primary  equations  used 

to  obtain  the  data  required  for  a  carpet  plot  design 

analysis.   Generally,  the  variables   V™  ,  DL ,  CT n    and 

T         LRo 

PT  ,  that  are  required  for  solution  have  specific  ranges 
of  values,  depending  on  the  weight  class  of  the  helicopter 
being  designed.   The  graphical  results  of  equation  (3.11) 
for  tip  speeds  of  600  to  700  fps  and  mean  lift  coefficients 
between  .3  and  .7  are  illustrated  in  Figure  5.1. 

Both  the  Fortran  and  Disspla  programs,  as  well  as  a 
decision  making  flow  chart  are  provided  in  Appendix  C 
to  aid  in  using  this  method  for  a  design  solution. 

E.   WEIGHT  EQUATIONS 

Weight  equations  need  to  be  developed  that  realistically 
reflect  the  sizing  class  of  the  helicopter  being  designed. 
The  evolution  is  greatly  simplified  if  a   specific  engine 
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installation  [#  and  horsepower]  is  assumed,  since  the 
weight  of  a  number  of  components  depend  only  on  the 
installed  power;  this  would  include  such  terms  as  the  engine 
controls  and  accessories.   Another  category  would  be  those 
components  whose  weights  depend  on  either  the  gross  weight 
on  two  or  more  of  the  following  in  combination:  rotor  tip 
speed  (VT) ,  rotor  diameter  (R) ,  rotor  solidity  (a) . 

The  equations  developed  here  are  taken  from  the  Hiller 
Aircraft  Corporation  Performance  Data  Report.  [Ref.  5]   In 
this  report  they  assumed  a  specific  engine  installation, 
the  Allison  T-63  with  a  military  power  rating  at  sea  level 
of  250  horsepower. 

There  is  a  possible  problem  of  the  validity  of  these 
weight  relationships  when  applied  to  different  helicopter 
design  categories.   However,  assuming  a  specific  engine 
determines  a  number  of  the  component  weights,  and  thus 
minimizes  the  inaccuracies.   Using  the  weight  estimation 
relationships  developed  in  the  Helicopter  Design  Manual 
[Ref.  2],  the  engine,  control  and  accessory  weight  can  be 
calculated  and  the  weight  formulas  developed  here  applied 
to  give  a  representative  useful  load  and  empty  weight 
formula  for  preliminary  design  analysis.   This  is  done  in 
Appendix  C,   for  a  20,000  pound  class  helicopter. 
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The  following  relations  are  used  to  reduce  the  component 
weight  formulas  for  the  specification  helicopter: 


W/DL  =  A  =  ttR' 


(3.12) 


W/PL  =  MHP  =  2  50 


(3.15) 


(Military  rating  for  Allison  T-63  at  sea  level.)   (PL  = 
Power  Loading.) 


P  =  /A7V, 


(5.14) 


Using  these  equations  the  component  weight  for  the 
specified  helicopter  empty  weight  may  be  reduced  to  the 


following 


Engine,  Controls  and  Accessories 


617.5  lbs 


Engine  Section 
Group 


053 [W/PL] 


1.07 


19.5  lbs.  (5.15) 


Main  Trans- 
mission 


10.43 


T(,1.29  5 

IV 


(PL  VT) 


863 


80 


1221  p-0UJ   (5.16) 


Rotor  Drive 
Shaft 


Tail  Rotor 


5.56 


W 


1.05 


(PL  VT) *  (DL) 


55 


52.22 


W 


1.14 


(PL  VT) 


1.7 


266  p-7     (5.17) 


17449 


V, 


1.14 


(5.18) 
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The  engine,  controls  and  accessories  category  includes 
such  items  as  lubrication  and  oil  cooling  system,  engines, 
communications,  engine  controls,  engine  accessories, 
instruments  starting  system,  furnishing,  flight  controls, 
electrical  system  and  stabilization.   These  are  considered 
fixed  weight  items  determined  from  specification  of  the 
engine  and  weight  class  of  the  helicopter. 


Tail  Rotor  W75 

Gear  Box    3.7  S «-  =  59.47  /P  (3.19) 

(PL  V T   (DL)  "" 


Tail  Rotor  w1.355  q7 

Drive  .124  2-— y-^  =  2.886  P     /A     (3.20) 

Shaft  (PL  VT)    (DL) " 

Body  and  916  qq 

Gear  =  1.91  W     +  .0294  W 
Landing 


Rotor  ..1.185    .33  ,  q ,-       -- 

Blade  j>5  . 1 5   t-^-f   =    o5.15   tt-  A  a  (.3.Z1J 

Teetering  V-CDL)*103  VT 


Rotor   Blade  1.205    .33  9nQ      _- 

Artie-  19.77    rrrr   =    19.77   rr-  A  a  (.3.ZZJ 

ulated  VT(DL)^UD  VT 


Rotor  Hub        wl . 21  ?1 

Teetering   .0088 ^r  =  .0088  WA  (3.23) 

DL'Zi 


Rotor  Hub         w1-21  71 

Artie-      .00975 j±    =  00975  WA*  (3.24) 

ulated  DL* 
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Fuel  System     .416  per  gallon  capacity  =  .0615  Wn     (3.25) 

r 

where   W_  =  fuel  weight. 


The  individual  component  weights  may  now  be  combined 
into  a  single  expression  for  the  helicopter  empty  weight 


W   =  617.5  +  .0617W„  =  1221P'863  +  266P*7  + 


F        x        "™        1.14 

T 

1     (3.26) 

+  58.47/P  +  2.886P,57/A  +  .191W916  +  .0294W'" 


+  appropriate  rotor  blade  and  hub  weights. 

As  stated  earlier,  the  design  specifications  called  for 
a  useful  load  consisting  of  a  pilot  (200  lbs),  payload 
(400  lbs)  and  the  required  fuel  weight  (W_) .   The  fuel 

r 

weight  is  calculated  for  the  Allison  T-63  in  the  following 
manner:  endurance  of  three  hours  at  85  percent  of  normal 
rated  powered  for  the  T-63  is  180.2  HP  and  the  specific 
fuel  consumption  at  this  power  is  .783  lbs  fuel/BHP  HR. 
Including  an  allowance  for  a  three-minute  warm-up  at  NRP 
and  using  a  5  percent  correction  factor  on  SFC,  as 
specified  in  Reference  5,  the  fuel  weight  becomes: 

Wp  =  3(180.2) (.822)  +  |g.  (212) (777)       (3.27) 

An  allowance  should  also  be  made  for  oil  plus  trapped 
fuel.   This  is  estimated  at  20  lbs. 
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The  total  useful  load  is  the  sum  of  the  useful  load 
items . 

Wu  =  200  +  400  +  452.6  +  20  =  1072.6  lbs     (3.28) 

A  new  variable,   WgAR  ,  is  defined  as  the  sum  of  the 
emply  weight  plus  useful  load.   It  is  the  of  equations 
(  3-  26)  and  (  3.  28)  ■ 


WBAR  =  1717.9  +  1221P'863  =  266P'7  +  17^4^4  +  58.47/P 

VT 


(3.29) 


57  —  916  99 

+  2.886P    /A  +  .191W     +  .0294W 


+  appropriate  rotor  blade  and  hub  weights. 

Equation  (3.11)  together  with  equation  (3.29)  form  the 
basis  of  a  carpet  plot  design  study.   These  equations  are 
solved  simultaneously  for  WBAR  •   This  solution  is  best 
illustrated  graphically,  as  in  Figure  5.2.   The  graph  in 
Figure  3.2  was  generated  for  a  specific  value  of   C. R   over 
a  range  of  tip  speeds- [600  to  700]. 

F.   GRAPHICAL  ANALYSIS 

Graphs  similar  to  Figure  3.1  are  generated  for  several 
value  of   CyR  ,  and  are  then  cross  plotted  to  form  Figure 
3.2. 

The  mean  lift  coefficient,  CyR  ,  values  are  selected 
based  on  what  is  considered  the  historical  average  range  of 
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values.   Figure  3.5  is  basic  plot  for  a  carpet  plot  design 
study.   Programs  are  provided  in  Appendix  D  which  will 
generate  the  required  data  sets  and  plots  of  Figures  5.2 
and  5.5. 

The  solution  field  depicted  in  Figure  5.5  is  too  large 
to  be  of  great  analytic  value  and  as  such  must  be  reduced. 
Three  parameters,  maximum  gross  weight,  rotor  diameter 
(both  specified  in  the  Design  Specification)  and  the 
aspect  ratio  can  be  used  to  narrow  the  field  of  solutions. 

1 .   Rotor  Diameter  Boundary 

A  net  to  exceed  value  for  the  rotor  diameter  is 
generally  given  in  the  design  specifications.   This 
limiting  value  is  based  on  the  operating  environment  of  the 
helicopter.   IVith   R  max  specified,  there  is  a  linear 
relationship  between  the  disk  loading  and  the  gross  weight. 


W     W 
DL  =  ?■  = 


A     D2 
t:R 


The  resulting  bracketing  of  the' solution  field  by 
applying  both  the  maximum  gross  weight  and  maximum  rotor 
diameter  limits  to  the  carpet  plot  are  shown  in  Figure  5.4 

2 .   Respect  Ratio  Boundary 

It  is  evident  that  a  further  restriction  is  still 
necessary  to  completely  define  the  region  of  acceptable 
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design  solutions.   Studies  have  indicated  that  a  main 
rotor  aspect  ratio  of  21,   is  a  representative  upper  limit 
Thus 


21  > 


R<mr> 
C<mr> 


b_ 

TTCJ 


b  po  CLR  VT2 

QTT   DL 


or 


DL  > 


b  p  CLR  V 
126tt 


For  the  case  of  a  two  bladed  main  rotor  equation  (3.50) 
reduces  to: 


DL  >  .000012  CLR  VT" 


The  detemination  of  this  boundary  graphically  is  as 

follows : 

2 
The  hover  solution  plot  of  Figure  5.2  is  replotted 

relative  to  the  coordinates  disk  loading  and  design  mean 

blade  lift  coefficient.   The  limiting  curves  for 

DL  =  .000012  CLR  \'2      are  then  plotted.   The  intersection 

with  the  appropriate  constant  tip  speed  lines  of  the  hover 

solution  represent  the  aspect  ratio  boundary;  Figure  5.5. 


For  a  helicopter  rotor,  the  aspect  ratio  is  defined  as 
the  radius  divided  by  the  chord. 

2 
For  clarity  lines  of  constant  gross  weight  are  omitted 
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These  intersection  points  are  then  cross  plotted  onto 
Figure  3.4.   Figure  3.6  represents  a  graphical  plot  of  the 
solution  set  satisfying  the  performance  and  structural 
design  criteria  of  a  small  observation  helicopter  as 
specified  in  this  study. 
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IV.   HESCOMP 

A.   DESCRIPTION  OF  PROGRAM  . 

HESCOMP  is  a  helicopter  sizing  and  performance  computer 
program  developed  by  the  Boeing  Vertol  Company.   The  program 
was  originally  formulated  to  provide  for  rapid  configuration 
design  studies. 

A  number  of  programming  options  are  available  to  the 
user  of  HESCOMP.   When  the  type  and  mission  profile  of  the 
helicopter  are  known,  HESCOMP  may  be  used  to  size  the 
aircraft.   Alternately,  it  may  be  used  for  mission  profile 
calculations  when  the  sizing  details  [gross  weight,  payload, 
engine  size,  etc.]  are  specified.   A  combination  of  these 
two  options  is  also  available;  the  program  may  be  used  to 
first  size  a  helicopter  for  a  primary  mission  and  then 
calculate  the  off-design  performance  for  other  missions. 
Finally,  HESCOMP  may  be  used  solely  for  obtaining  helicopter 
weight . 

Sensitivity  studies  involving  both  design  and  per- 
formance tradeoffs  can  easily  be  done  with  HESCOMP.   Incre- 
mental multiplicative  and  additive  factors  can  be  imbedded 
in  the  input  data. 

The  various  helicopter  configurations  that  may  be 
studied  using  HESCOMP  are  detailed  in  Table  4.1. 


54 


TABLE    4.1 


HELICOPTER  CONFIGURATIONS 
WHICH   MAY    BE    STUDIED    USING   HESCOMP 


HELICOPTER  CONFIGURATIONS  WHICH  HA*  BE  STUDIED  USIHG  HESCOMP. 

^~"~--^Addltlonal  Lif  t/Propul.i  ion 
'  -~-S.y3tcm  Components  Which 
~^~"- — JJunt  be  Added  to 
Helicopter           ■ — ^"Pure" 

Type                  ^-\Conf. 
(Doth  Single  I  Tandem  Rotor)    ^~"-^_^ 

Wing 

Propeller 

for 
Auxil  lary 
Propul slon 

Auxi llary 

Independent 

Engines 

Type  of  Auxil lory 
Independent  Engines 

1 

T/Shaft 

T/Fan 

T/Jet 

Pure  Helicopter 

. 

winged  Helicopter 

X 

Compound  Helicopter 

(1)  Coupled  (prim,  engines 
drive  auxiliary  propulsion 
sy9 tern) 

(2)  Auxiliary  Independent 
propulsion  system 

(u)  T/Shaft  engine 
(b)  T/Fan  engine 
(e)  T/Jet  engine 

X 

X 

X 
X 
X 

X 

X 
X 
X 

X 

X 

X 

Auxiliary  Propulsion  Helicopter 

(1)  Coupled  (prim,  engines 
drive  auxiliary  propulsion 
system) 

(2)  Auxiliary  Independent 
propulsion  system 

(a)  T/Shaft  engine 
(b|  T/Fan  engine 
(c)  T/Jet  enalna 

X 

X 

X 
X 
X 

X 

X 

X 

oaxial  Rotor  Helicopter 

(1)  Coupled  (prim,  engines 
drive  auxiliary  propulsion 
system) 

(2)  Auxiliary  Independent 
propulsion  system 

(a)  T/Sliaft  engine 

(b)  T/Fan  engine 

(c)  T/Jot  engine 

X 
X 

X 
X 

X 

X 

X 

X 
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B.   PROGRAM  MODIFICATIONS  AND  IMPLEMENTATION 

The  computer  program  received  from  Boeing  Vertol 
required  some  modification  and  reformating  in  order  to  run 
properly  on  the  Naval  Postgraduate  School  IBM  system. 
These  alterations  did  not,  however,  alter  the  program  output 
or  usability. 

HESCOMP,  as  received  from  Boeing  Vertol,  was  17821 
lines  long  and  set-up  as  a  sequential  data  set  to  be 
assemble  on  a  'G  compiler'.   The  Batch  processing  system  at 
the  Naval  Postgraduate  School  accepts  only  programs  set  to 
run  on  'H  compiler'.   Normally,  the  differences  between 
these  two  compilers  are  minor  and  programs  that  run  on  one 
will  run  on  the  other.   However,  this  was  not  the  case  with 
HESCOMP. 

In  order  to  facilitate  the  program  debugging  process, 
HESCOMP  was  reformatted  as  a  partitioned  data  set.   What 
this  effectively  did  was  to  break  the  program  down  into 
eight  members  of  approximately  2000  lines.   The  program 
breakdown  is  illustrated  in  Table  4.2. 

Each  of  these  were  compiled  individually  and  then  error 
codes  analyzed.   The  member  data  set  was  then  modified  as 
required  to  properly  compile. 

Once  all  the  members  of  the  partitioned  data  set 
compiled  properly,  HESCOMP  was  again  formated  as  a 
sequential  data  set  and  run  utilizing  input  data  for 
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which  there  was  a  known  output.   This  insured  that  the 

modifications  made  to  the  original  program  had  not  altered 
the  logic,  ie.,  gave  faulty  results. 

The  control  language  program  to  access  HESCOMB  on  the 

Batch  processing  system  and  a  sample  input  and  out  data 
set  are  shown  in  Appendix  D.     These  are  also  available 
on  the  Aero  disk  for  copying  and  use. 

TABLE  4.  2 

PARTITIONED  DATA  SET 


MEMBER  NAME 


LINE  NUMBER 


SIZE 


FIRST  ROUTINE 


SI 
S2 
S5 
S4 
S5 
S6 
S7 
S8 


1  - 

■  1681 

1681 

AERO 

1682  - 

•  4132 

2451 

CLIMB 

4153  • 

•  6551 

2399 

XIBIV 

6532  - 

■  8974 

2443 

POWAVL 

8975  ■ 

•  10870 

1896 

PRINT  1 

10871  ■ 

■  13042 

2172 

ROT  POW 

13043  ■ 

-  15583 

2541 

CRUS  5 

15384  ■ 

■  17821 

2448 

TAXI 

C.   PROGRAM  FLOW 

The  program  is  conceptually  outlined  in  Figure  4.1, 
[Ref.  7].   The  program  flow  is  monitored  by  a  general  loop 
which  controls  a  series  of  peripheral  programs.   There  are 
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INUNl    AMU  PHOflll  EH/KUIOR    QATA 
.  NLINI     SIZING    LSI  It  MIA 
VllMUl    UOHIRIC    Of  f  INI  T  IONS 
SUH'.TSIlN   tf  ICliT    1RINI)   COIFFICIFNTS 
VI  Ml  CI  L    AtRUUTMNIC    UIASAUIRISIICS 
MISSION    PHUl  lit    OAIA 


MISSION 

PERFORMANCE 

SUBPROGRAM 


ACCESSES    9 

PERFORMANCE 

SUBROUTINES 


INPUT 


AERO 

OTNANICS 

SUUHOUflNE 


SWE    1R1NOS 
SUOROUIINt 


iNGINE 
SUING 
SUURDUlINf 


Mi  IGKIS 
IRtNOS 
Su« ROUTINE 


OOII 


U'^ 


1. 

OflAlltO  VtHlCll     OININMON'j 

1  • 

SUUSTSKMS   Mt  IGHI    aWAADU* 

1  t 

ENGINE    BAKU    PUWI  R 

1    • 

PARASIIt    DRAG    BRIArDOUN 

NISSION    PIRtOHMANU     UNI    H»TOnT 

Figure    4.1.      HESCOMP    Program    Flow 
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a  total  of  44  subroutines.   Detailed  program  descriptions 
cam  be  found  in  Section  4  of  the  HESCOMP  User's  Manual. 

D.  PROGRAM  INPUT 

Program  input  can  be  loosely  group  into  ten  categories : 
general  information,  aircraft  descriptive  information, 
mission  profile  information,  rotor  tip  speed  schedule, 
incremental  rotor  performance,  auxiliary  propulsion  input 
schedule,  engine  cycle  information,  rotor  performance 
information,  propeller  performance  information,  and 
supplementary  input  information. 

The  actual  amount  of  input  data  requires  varies  greatly 
with  the  program  options  selected.   An  example  of  a  data 
set  formatted  to  run  on  the  IBM  system  is  shown  in 
Appendix  E.   A  more  detailed  explantion  is  available  in 
Section  5  of  the  HESCOMP  User's  Manual. 

E.  PROGRAM  OUTPUT 

An  example  of  the  program  output  is  included  in 
Appendix  E.   The  printout  consists  of  general  data,  input 
data,  sizing  data  [program  output]  and  mission  performance 
data  [for  the  size  helicopter] .   Detailed  descriptions  of 
these  and  diagnostic  error  statements  are  described  in 
Section  6  of  Reference  6. 


59 


V.   CONCLUSIONS  AND  RECOMMENDATIONS 

Three  approaches  to  analyzing  a  preliminary  helicopter 
design  were  explored  in  the  course  of  this  paper.   It  was 
found  that  a  number  of  the  performance  equations  could  be 
greatly  simplified  with  little  degradation  in  the  final 
results.   A  sensitivity  analysis  brought  further  insight 
into  the  inter-play  of  the  parameters  and  how  changes  in 
them  tended  to  effect  the  helicopter  performance  equations. 

Carpet  Plots  provided  the  most  interesting  method  of 
analysis.   Development  of  a  graphical  solution  matrix 
using  this  method  provides  a  usual  interpretation  of  what 
is  occurring  when  key  parameters  are  varied. 

Two  cases  were  explored;  a  light  observation  helicopter 
in  the  5,000   pound  weight  class  and  a  heavier  utility 
helicopter  in  the  20,000  pound  weight  class.   The  Carpet 
Plot  method  provided  reasonable  solutions  in  both  cases. 
In  doing  the  analysis  for  the  utility  helicopter,  the 
initial  weight  estimation  equation  had  to  be  adjusted 
upward  by  approximately  2,000  pounds  for  the  equations  to 
intersect  properly.   This  is  not  considered  a  limitation 
to  this  method  of  analysis,  however,  it  does  point  up  an 
area  for  further  investigation.   It  may  be  possible  to 
develop  more  accurate  weighing  factors  for  this  equation 
when  dealing  with  higher  gross  weight  helicopters. 
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HESCOMP  provides  a  plethora  of  information  to  the  user 
However,  the  price  is  the  amount  of  inputed  data  required 
for  even  a  simplified  analysis.   At  a  preliminary  design 
level  of  analysis,  the  other  methods  explored  provide  a 
quicker  first-cut  look  at  the  potential  design. 
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APPENDIX  A:   NOMENCLATURE 


TERM  DEFINITION  UNITS 

a  Slope  of  Airfoil  Section  Radians 

Lift  Curve 

2 

A  Rotor  Disk  Area  ft 

AR  Aspect  Ratio  Dimensionless 

ATR  Tail  Rotor  Disk  Area  ft2 

b  Number  of  Rotor  Baldes  Dimensionless 

B  Tip  Loss  Factor  Dimensionless 

C  Main  Rotor  Cord  ft 

C,  Profile  Drag  Coefficient  Dimensionless 

at  Zero  Lift 

C,  „  Design  Mean  Blade  Lift  Dimensionless 

Coefficient  at  Sea  Level 

CT  Coefficient  of  Thrust  Dimensionless 

C  Coefficient  of  Power  Dimensionless 

P 

5  Blade  Section  Drag  Dimensionless 

Coefficient 

DL  Disk  Loading  lb/ft2 

FM  Figure  of  Merit  Dimensionless 

HP  Horsepower 

LTR  Tail  Rotor  Moment  Arm  ft 

2    4 

p  Air  Density  lb  sec"/ft 

y  Advance  Ratio  Dimensionless 

R  Rotor  Radius  ft 
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TERM  DEFINITION  UNITS 

PT  Total  Power  HP 

PTM  Main  Rotor  Total  Power       HP 

PTtd  Tail  Rotor  Total  Power       Hp 

P^  Profile  Power  HP 


P-  Induced  Power  HP 

P  Parasite  Power  HP 
P 

PL  Power  Loading  LB/HP 

R  Rotor  Radius  ft 

T  Thrust  HP 

V,  Induced  Velocity  ft/sec 

Vp  Forward  Velocity  ft/sec 

V  Aircraft  Forward  Speed  ft/sec 

VT  Rotor  Tip  Speed  ft/sec 

W  Aircraft  Gross  Weight  lbs 

W  Empty  Weight  lbs 

W-  Fuel  Weight  lbs 
F                    5 

W  Useful  Load  lbs 
u 

WBAD  Empty  Weight  Plust 

BAK  Useful  Load  lbs 


Solidity  Dimensionless 
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APPENDIX  B:   CARPET  PLOT  FORMULATION  FOR  2  0,0  00  LB 
CLASS  HELICOPTER 


Bl     SPECIFICATIONS: 

Maximum  Gross  Weight:   20,000  pounds 
Maximum  Rotor  Diameter:   30  feet 

B2     PRELIMINARY  ENGINE  SIZING: 

B2 . 1   Utilize  equation  (2.14)  to  determine  engine 
horsepower  category. 


W  =  [4.753PTR]2/o 


20,000  =  [47.53PT  30]2/3 


PT  =  1985  HP 

B2.2   Use  the  engine  selection  parameters  tables  B.l 
to  determine  the  number  and  type  of  power  plant 
[table  taken  from  Reference  3] . 

B2.2a   Type  and  number  selected:   2  type  C. 

B2.2b   Specifications: 

Dry  Weight  Per  Engine:   423  pounds 

Shaft  Horsepower  at  Standard  Sea  Level: 

Military       1561  HP 

Normal        1318  HP 

B3     WEIGHT  EQUATION  FORMULATION 

B3.1   To  obtain  the  engine  control  and  accessory 
weight  use  items  7,  9,  10,  11,  12  and  13  of 
the  weight  estimation  relationships  developed 
in  Reference  5  for  a  utility  helicopter: 
#7:  609  lbs;  #9:  129  lbs;  #10:  76  lbs; 
#11:  410  lbs;  #12:  459  lbs;  and  #15:  502  lbs. 
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TABLE  B.l 
ENGINE  SELECTION  PARAMETERS 

The  following  turboshat  power  plant  data  are 
presented  for  one  engine. 

Engines:  A     B     C     D*    E     F 

Dry  Weight  (lbs)     158    288    423    709    580    750 

SHP  (ssl)  Military  420  708  1561  1800  2500  3400 
Normal  370  659  1318  1530  2200  3000 
Cruise     278    494   1989   1148   1650   2250 

SFC  (ssl)  Military  .650  .573  .460  .595  .615  .545 
Normal  .651  .573  .470  .606  .622  .562 
Cruise    .709   .599   .510   .661   .678   .610 

Initial  Costs       $93K   $100K  $580K  $560K  $640K  $700K 

Operating  Cost 

per  hour/engine      $8     $16    $20    $35    $40    $60 


Preventative  Maint 


per  hour/eng 

ine 

$25 

$50 

$1 

MTBMA  (hrs) 

3.5 

3.0 

2.0 

MDT  (hrs) 

0.7 

0.6 

0.5 

MTBF  (hrs) 

185 

210 

205 

MTBR  (hrs) 

600 

750 

800 

$50   $100   $125   $160   $220 

3.0    4.0  3.5 

1.3    2.0  2.6 

285    280  320 

800   1000  750 
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B3.2    Simplifications 


W     ,     D2 

DL  £pm 


W  'A 

=  MHP  =  31,00  ;  p    =      ±- 


B5.5   Engine  Group 


.055(5100)1,07  =  272  lbs 


B3.4    Main  Transmission 


10.43 


IV 


1.295 


Upm   Vt) 


163      ~    .432 
LAJ 


=    10.43 


W865   A+-432 

Upmr863vT-863 


(10.43)  (3100)  •86-5p-863 


=    10,748P 


863 


B3.5         Rotor    Drive    Shaft 


5.56 


W 


1.05 


.7      w    .35 
(£pm   VT)  [£] 


=    5.56(5100)  'V7 


=    1545P 


.7 


B3.6        Tail    Rotor 


32.22 


W 


1.14 


(Jlpm   VT) 


1.14 


307,600 

V   T7TA 

VT 
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Tail  Rotor  Gear  Box 


.  75 


3.  7 


= ~   =  (3.7)(3100)   P 


o  ^  .   o 


w 


(4  pm  VT)    [£] 


=  206P 


.5 


3.8    Tail  Rotor  Drive  Shaft 


124 


W 


1.355 


Upm) 


5  7  .  r . 7  8  5 

w 


(.124)  (3100)  -57P'57  /a 


A 


85 


12.12P'57  /A 


3.9    Landing  Gear 


916  99 

=  .191Wyi°  +  .0294W 


3.10   Rotor  Blades  Articulated 


19.77 


,.1.206  .33 

W     a 


VT  DL 


.205 


.205  .33 
=  19.77  ~  A  a 

vT 


B3.ll   Rotor  Hub  Articulated 


00975 


IV 


1.21 


00975WA 


.21 


DL 
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B5.12   Fuel  System 


0615  W, 


Calculation  of  fuel  weight  three  hours  at 
cruise  SHP 
1513  lbs  +  10% 
1664  lbs 

3.13   Total  Equation 


WB  =  12,987,*  +  107948P-863  +  1545P"7 


307600    in,n.s  c7 

+    J  14  +  206P,:>  +  12.12p-b7  /A" 

vTi-14 


916  99 

+  .191W      +  .0294W 


W     9  0  5   3  3  21 

+  19.77  ^-  A'"u:5S     +  .00975WA 

VT 


B4    HOVER  EQUATION 

Following  the  formulation  in  Section  of  Chapter  3, 
the  weight  equation  based  on  the  design  mean  lift  coeffi 
cient  and  power  required  is: 


K, 


W  = 


1  ■ 

■    411 

51 

DL3/4  I 

1 

+ 

K7 

vT 
/ur 

1/2 


-  K 


VT  +  K5  /DL 


This  number  was  increased  from  8987  to  12987  to  bring 
the  curves  together.   This  reflects  a  4000  lb  useful  load. 
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where : 


9583  ? 

=    -^^    (i    +    1.8073    CTD    n 

1  CLRo  LRo 


K2         PT6000/950    1J£-1 

Kl 


0.00025929       ,,         ,     ,„„,.„         2. 


K 

'"LR 


C~ CI    *    l-8078CLRb-) 


o 


„  553480.0 

K4  =  — r 


v  5695.7 

K5    =         K1 


B. 5   GRAPHICAL  RESULTS 

Figure  B.l  is  an  example  of  equation  (3.13)  plotted 
against  equation  (B.4)  for  a  specific  design  mean  lift 
coefficient . 

Figure  B.2  illustrates  the  family  of  curves  obtained 
when  the  design  mean  lift  coefficient  is  varied  from 

0.3  to  0.7  . 

In  Figure  B.3  the  solution  matrix  depicted  in  Figure 
B.2  is  narrowed  by  the  constraints  placed  on  the  gross 
weight,  rotor  diameter  and  aspect  ratio. 
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APPENDIX  C.   CARPET  PLOT  METHODOLOGY  FLOW  CHART  AND 
EXAMPLE  PROGRAMS : 

This  section  contains  a  flow  chart  to  help  organize 

a  carpet  analysis  and  example  IBM  computer  programs  to 

produce  the  data  sets  and  disspla  graphs. 
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C ******** *************44444********<4*****44444************************* 

»0, ,.,»,».               GRAPHICAL    fcELICCPTEF    DESIGN    EFOGRAM  ********** 

c********                                     CARErT    PLCT     NUCEER    1  ******** 

«,,,...»                                             EY    AL     HANSEN  ******** 

C********  ******** 

c********            TEIS    pscGFA*     IS    DESIGNII    TO    PBOVIDE    A    HEIGHT    VS.  **4***»* 

........                  CISK    ICAEING    GRAPH     Al    SEVERAL    DIFFERENT    ROTOH  ******** 

........                  TIP    VELOCITIES     FOR     A    SPECIFIED    DESIGN    SEAN  ******** 

C***»****                                                LIFT    COEFf  fICEN'T.  ******** 

C**********                                             OBSERVATION    CLASS  ********** 

C***4**************** 444** 4********4****** ****************************** 
C******»*************444444***************************************4»**4* 

c****  ***** 

C**                                                                         NOKENCIATURE  *** 

C*  ** 

C*     VJBIABLES:  ** 

C*  ** 

C*     CLE      EESIGN  KEAN  LIFT  COEFIICENT  ** 

C*     VI       TIP  VELOCITY  ** 
C*     DL       DISK  LCAEING 

C*     b         WEIGHT  IS    CALCULATED  FFC."i  POWER  EQUATION  ** 

C*     i'E        CESFOL  LCAD  ILUS  EMPTY  WEIGHT  ** 

C*     Pi       FCWEE  AVAILAELE  IN  HORSEPOWER  ** 

C*  *» 

C****  4**4* 

£4*44*********************4****4***44*********************************** 

31 

1  KJ 


►  *******************4****4***<  4*****  **»******************: 

3EAL    CLR,EL  f  1  5C )  ,W1  (150)  .  W2  (1  5C)  ,  W3  (150)  .  W4  (150)  ,  W5  (150)  , 
•  E1  (15  0)  ,WB2  (  150)  ,WE3  (IbO)  ,  WEM1:>0)  ,  «E5  (150) 


INTEGER    I 
C 

c 

C DEFINE    FILES 


C 


CALL    FBTCMS  ('FILiDEE    ','3  ','DISK             •/,CR?T35' 

1  'EATA             ■  ,  •  A  ') 

CALL    FRTCHS  ('JILEEEr    ','U  ','CISK             ','CRPT45' 

1  TATA             •      •  A 


C READ    DATA    FROM    FILE:    C2PTX    DATA    A    

C 

EC  10  1=1,99 
C      CRPT1  DATA  A 

HEAD  (3,"/C)  DL  (I)  ,Vi  1  (I)  ,  WE1  (I)  ,W2  (I)  ,WE2  (I) 
C      CRPI2  DATA  A 

BEAD  (4,71)  W3  (I)  ,WE  2  (I)  ,  W4  (I)  ,WB4  (I)  ,  W5  (I)  ,  WB5  (I) 
10     CONTINUE 
C 

C CALL  DISSPLA  RCCTINES  FOB  PLOT  

C 

CALL  TEK618 

CALL  MECEUF 


CALL  RESET  (3H2IL) 
CALL  HWSCAL  ('SCREEN') 
CALL  PAGE  (12-C,  9.5) 


:l.i 


,.H, 


CALL  GBACE  (0 

CALL  PHYSCR  ( 

CALL  AEEA2D  (10.0,6. 

CALL  FBAKE 

CALL  SklESL 

CALL  BASALF  ('L/CSTC') 

CALL  MIXAIF  (  'STAND') 

CALL  INTAXS 

CALL  YTICKS  (5, 

CALL  XTICKS  (5 


CALL  SHDCHR  (  -  90 ,  1  , . 015 , 1 ) 
CALL  HEIGET  (-  1b) 


CALL  XNARE"  (''  (D)  ISK  (L)  0  A  DI  NG  2  ' ,  1  00) 

(W)  EIGHT     |  (L3S)   ) 


CALL    YNAKE     ('(G)RCSS     (W)  EIGHT     (  (L3S)  )  $  •  ,  1  00) 


CALL    HEIGHT  (.290) 

CALL  HEAEIN  ( '  J  H  )  ELI  COPT  E  R  (C)i>RFET   (P)LCTS:   (CLR=0.5)S' 

C      CALL  MESSAG  ( '  < H )'  ELI  cCPTE  5  (C)£EF£I  (P)LCTS:  (CLR=0.5)3' 

C     1  100,3.25,7.55) 


75 


70 
71 


CALL    H 

CALL    T 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL    B 

MAXLIN 

CALL    H 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

STOP 

FORMAT 

FOEMAT 

FOEHAT 

END 


EIGET     ( 

HAJ-       il. 

HKCEV 
EGLIN 
URVE 
URVE 
URVE 
ORVE 
DEVI 
ASE 
URVE  (D 
U8VE  (D 
URVE  (D 
URVE  (D 
URVE  (D 
HKCEV  ( 
LREC  (4 
HIE  (1, 
LOEF  (1 
=  LIKEST 
EIGHT  ( 
INESP  ( 
IMS  (' 
INES  (' 
INES  (• 
INES  (' 
YLEGN  ( 
EGEND  ( 
ONEEL 


.  20) 

j.-  i,2.9,225C.  ,5<j.  .zoou 

(.C18) 


(DI  ,S1  ,95  ,0 
[ULf  S  2,95,0 

B" 

(D 


I,»3,95,0) 
L,H4,95,0J 

L,*S,95,0) 


L,H 
I  #8 

L,H 

.C3 
-2, 
1YG 
) 


E1,95,0) 

E2,95,0 

E3,95,0> 

54,95,0) 

=5,95,0) 

0) 

4.65,1.6  ,  1.  35,1) 

BID) 


1$«  ,IEAK1 


(IEAK1,300,20) 

2.0[ 

( '» )  E I G  H 

(I) NEUC 

(E FCFI 
<E) ABAS 
'  (?)  CHE 

IEAK1,4 


EDJ{  ,  IPAK  ^,2) 
LES'  , IPAK  1,3) 
II£J»  -IPAK  1,4) 
E     (C)  URVES1  ,16) 
,4.43  ,4.8) 


STATES 

(5  (2X, 

6(2X, 


EST! 

F  10. 
F  10. 


HI 


76 


c** 
c** 
c** 
c** 

c** 

c*» 

c** 

c** 

c** 

c** 

c** 

c** 

c** 

c*» 

c»* 

c** 

c** 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c** 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c** 


»**********»******4*4*******»***44*************************** 

********  GHAEHICAI    EELICCPTEE    DESIGN    E50GRAS  *  * 

******  FAK1LY    OF     SOLUTIONS  * 

******  CAEEET    PLOT    NUBEER    2 


****** 
****** 
****** 
****** 
****** 
****** 
****** 
****** 


CAEEET    PLOT    NUBEER    2 
ES    AL     HANSEN 


FAMILY 


C3SE2VATICN    CLASS 
* *******  ** 

4  ***  ****************  ************  *  ******  ********************** 
»4*********»********»***********4**************************** 

** 

NOMENCIATUEE 


»  *  * 

*  *  * 

*  *  * 

*  *  * 
**  * 
*** 

*  *  * 

*  ** 

*  *  » 
**  * 

*  *  * 
*** 

*  *  * 

*  *  * 

*  *  * 


VAHIA3LES: 

CLE  DESIGN    MEAN    LIFT    COEFEICENT 

VT  TIP    VELOCITY 

DL  CISK    LOADING 

i  WEIGHT     AS    CALCULATED     EECM    POWEE 

BE  CSEFUL    LCAD    E  LU  S     EMPTY    WEIGHT 

EA  ECi.EE     AVAILABLE    IN    HCBSEPOWEB 


EQUATION 


C 
C~- 


** 

DL3 
EL4 
DL5 
CL6 
CL7 
V3 

an 

V5 
■  0 
17 

BVT1 

HVT2 
WVI3 
WVT4 
SVT5 
EVT1 
DVI2 
CVT3 
DVTU 
£VI5 

4  ********  ** 

BEAL*4  D 
*W5(5)  ,i6 
*.VT5(5)  , 


***  ** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 

*  *  *  * 
***** 
***** 
***** 
***** 

*  *  *  ** 

*  ** 
** 
** 
** 
** 
** 

** 
** 
** 
** 
***** 


EQUA 

EQUA 

EQUA 

EQUA 

EQUA 

EQUAL 

EQUAL 

EQU  AL 

EQUAL 

EQUAL 

"     EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EQU 

EvU 


IS  TH 

LS  Tb 

LS  TH 

LS  It 


LS  TH 
S  THE 
S  THE 
S  THE 
S  THE 
S  THE 
ALS  « 
ALS  Vi 
ALS  W 
ALS  W 
ALS  W 
ALS  T 
ALS  T 
ALS  T 
ALS  T 
ALS  T 
***** 

L3<5) 
(5)  .5; 
CVTl  ( 


E  DISK 
E  CISK 
E  DISK 
E  CISK 
E  DISK 
WEIGHT 
WEIGHT 
SEIGHT 
WEIGHT 
WEIGHT 
EIGHTS 
EIGHTS 
EIGHTS 
EIGHTS 
EIGHTS 
EE  CCER 
HE  CCRS 
fcE  CCEE 
HE  CCRE 
EE  CCRE 

*  *  4  4  *  *  * 

ELu  (5 


LOADING     FCE 

LOADING    FCF 

LOADING     FCE 

LOADING     FCE 

LOADING    FCE 

FOP.    CLR=.2 

FOE    CLB=.4 

FOE    CLE  =.5 

CLS=.6 

CLB=.7 

600 


CLa= .3 
CLR=  .U 
CLR=  .5 
CLB=  .  o 
CLiv=  .  7 


FOR 

FOE 

at  v: 


.ELU  (5) 
7  (5)  .SV 

5) ,DVT2 


AT    VT=o25 

AT    VT=650 

AT    VT=675 

AT    VT=700 

ESPONDING  DISK    LCADING    AT    VT=600 

ZSPONDING  CISK    LCACING    AT    VT=625 

ESPON'DING  DISK    LCADING    AT    VT  =  650 

ESPONDING  CISK    LCACING    AT    7T=675 

ESPONDING  CISK    LCACING    AT     VT=700 

4********4  4*****  ********************** 

,DL5  (5)  fDLC(5)  ,DL7f5>  ,»J  (5)  ,*4(5)  , 
11(5)  ,HVT2{5)  ,WV13<5)  ,WV1<1  (a)  , 

(5)  ,  DVT3  (5) ,DVT4  (5)  ,  DVT5  (5) 


******** 


DEFIt. 

CATA 

DATA 

CATA 

LATA 

CATA 

DATA 

DATA 

DATA 

LATA 

DATA 

CATA 

DATA 

CATA 

DATA 

DATA 

DATA 

CATA 

DATA 

DATA 

DATA 


E  DATA  — 
DL3/2.06 , 2 
■3/2505- ,2 

DL4/2.42, 2 


Wo/2424. 66 
DL7/2.60.2 
S7/2411.92 
WV^  1/250  5. 
W VT2/2478. 
WVT3/2455. 
WVT4/2433. 
WVT5/241 S. 
DVT1/2.06  , 
DVT2/2.07, 
DVT3/2. 07  , 
DVT4/2.07, 
DVT5/2. 05  , 


6  .Z 

M,2 

44  2 
5,2 

4  19 

Ut5  2 

h2 
390 
467 
244 
,2* 
.23 

,  z  J 

«2, 
44, 
46, 

«7  , 
«7 


.C7 
,24$ 

.40, 

•95, 
.67, 
.51, 
.79; 

.43, 

.65, 
.72, 
.17; 
2-95 
19.o 
99.7 
62.9 
2.63 
z.65 
2.67 
2.69 
2.69 


2.07, 
5.46, 
2.47, 
2419. 
2.69, 
2  399. 
2.81, 
2383. 
2  .86, 
2  372. 
2  440. 
.24  19 
2,239 
2,237 
9,2  36 
,2.74 
;2.77 
,2.79 
,2.81 
,2.81 


,05 
I33 
,47 

.69 

,a1 

'■i'i 

.67 

i:3 

,13 
.65 

H 

:.  O 

;.8 
1.  a 

1.8 


35  ,2418.25/ 


399.72,2382.99/ 

379  .13,2304.65/ 

/ 

365.10,2352./ 

357.35,2342.25/ 
424.67,2411 .92/ 
2402.53.2390.72/ 
3*3  .58,2372.5/ 
,2365.  1.2357. 35/ 
,2352. ,2342.25/ 

% 

5/ 
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C CALL  DISSELA  RCCIINES  FOB  PLOT 

C 


c 

C70 
C71 


ALL 
ALL 
ALL 
ALL 
XLL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 


TEK618 
MEDEOF 
BESET     (3 
HSSCAL     ( 
PAGE     (12 
GFAC1     (0 
PHTSCB     ( 
AEEA2D     ( 
FBARE 
3WISSL 
3ASALF     ( 
dlXALF     ( 
INTAXS 
YTICKS     ( 
XTICKS     ( 
SHDCHB     ( 
HEIGET     ( 
XNABE     (' 
YNAKE     I' 
HEIGHT     ( 
HEAEIN     ( 


HALL) 

•  SCREEN') 

:S{  e-5» 

1.0.  1.2) 
10-0,6.5) 


'I/CSTE') 
•  STAl>E' ) 


50,1, .015,  1 


.  Ifai 

E)  I! 
(Gj  BCSS  '  (W)  EIGHT     ('(LBS)  )  $' ,  1 00) 

29f 


(ElISK     (L)  OAOINGS',  100) 


ALL 

ALL 

YGRID 

ALL    T 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL    G 

ALL    3 

AXLIN 

ALL    H 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

ALL 

TCP 

CHHAT 

OHflAl 

OBBAT 

ND 


EIGfiT 

B  A F     (2 

=  5 

HKCEV 

ABA3 

EGLIN 

URVE 

OBVI 

OHVE 

UHVE 

uhve 

ASfc 

ORVE 
DBVE 
DBVE 
UHVE 

uavE 

HKCBV 
LREC 
BIE  ( 
LOFF  (1 
=LINEST 
EIGHT  { 
INESP  ( 
INES  (' 
INES  (' 
INES  (' 
INES  (• 
YLEGN  ( 
EGEND  ( 
ON  EEL 


10) 

1  JH)  ELICCPTEfi     (C)fRPET      (P)LOTSS', 

ICC,  1.0,  1) 

.20) 

0,.  1,3.  C,230C.  ,50.  ,2550.) 


(.C18) 


[D 

)D 
( 0 
[D 
(D 

(1 

B 

( 


L3,W3,5,0 

ljj;bu;5;o 

L5,«5,5,0 
L6, 1(6,  3,0 
L7,»7,5,0 


V1 1  ,iVT1,5,0) 

VT2,»VT2,5, 0) 

VT3,.VT3,5,0) 

VT4,iVT4,5,  0) 

VT5.iVT5,5,0) 

.030) 

.2.^.65,1.6,  1.35,1) 

IYfiRID) 

) 


(IEAK  1  ,300,  20) 

.  12) 

2.0} 

(h)  EIGHTS'  ,IEAK1 

(I 

(E. 

h. 

IEAK1,4  ,4.4 


NEUCEDS1  ,  IP  A K  1,2 
HCFILE5'  , IPAK  1,3) 


AEA5ITE3  '  ,IPAK1,4) 

)  CWEB     (C)  UHVE  £'  ,16) 

3  ,' 


(IB  V  ES 
4  .8) 


STATEMENTS 
(5  {21,1  10.3 
6 (2X, F  10 


rill 
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BY    AL     HANSEN  *»*4***« 

**  ****** 
B  IS  DESIGNEE  TO  GRAPHICALLY  DETERMINE  ******** 
T     SATIO    BOUNEARY    REQUIREMENTS    FOR  ******** 

YSTEM    USING    E5EVIOUSLY    GENERATED    DATA  ******** 

********** 


~  *****  *  ************  *444444*****4***4  4*****  4444  4************************* 

c********«*  GBAtHICil    fcELLCCPTEE    DESIGN  PEOGRAM                               ********** 

;»«"""»  ASPECT    RATIC    50UNDARY  ***4»*** 

c********  L0CI    Cf    HCVER    AND    USEFUL-LCAD    SOLUTIONS  ******** 

c***4****  CAEEET    PLCT     NDCEER    3  ******** 

c* ,,,*,».  3Y    al     HANSEN  ***4**4* 

C***<****  ******** 

c**44***4  TBIS    PECGEAM 

«*.»»««.  THE     ASPECT     , 

c. **,.,.,  A    ECTCR    SYSTEM 

C**********  ***»,,»»»* 

c**44****************444*»*********4**4********************************* 

C**»******4t*********44444*4  ******  **44*********************************** 
C****  ***** 

C**  NOMENCLATURE  *♦* 

C*  ** 

C*  V8BIABLES:  ** 

C*  ** 

C*     CLE      DESIGN  MEAN  LIFT  COEFIICENT  ** 

C*     VT       TIP  VELOCITY  ** 

C*     EL       DISK  LCADING 

C*     AE       ASPECT  EATIC.  HISTORIC ALLLY  ASSUMED  TO  BE  LESS  THAN  21   ** 

C*  ** 

c****  ***** 

C*  DVI1    EQUALS    TEE    CCRR  5SPC  II  DI  NG    DISK    LCADING    AT    VT  =  600 

C*  CVT2    EQUALS    TfcE    CCFB ESPON EING    DISK    LCADING    AT    VT=625 

C*  DVI3    EQUALS    THE    CCRR ESPO N DI NG    DISK    LCADING    AT    VT=o50 

C*  EVT4    EQUALS    TfcE    CCRR ESPO N DI NG    DISK    LCADING    AT    VT=675 

C*  DM15    EQUALS    TEE    CORRESPONDING    EISX    LCADING    AT    VT=700 

C*  C600    EQUALS    TfcE    LIFT    COEF    AT    Vl=o00 

C*  C625    EQUALS    TEE    LIFT    COEF    AT    VI=o25 

C*  C650    EQUALS    TfcE    LIFT    COEF    AT     U=o50 

C*  C675    EQUALS    TfcE    LIFT    COEF    AT    VT=o75 

C*  C700    EQUALS    TfcE    LIFT    COEF    AT    VI=7G0 

C**»*****»*»***********************< ************************************ 

BEAL*U  CLE  (5)  , EL  (  10)  ,C600  (1  0)  ,C625(1  C)  ,Cb5u  (  10)  , 
*C675(1  0)  .C70J  (  10)  , 

*DVl1(5),D*I2(5),DVT3l5),DVT4(5)  ,DVT5  (5) 
C 
C DEFINE  DATA   

DATA  DVT1/2.  0b,2.  U2,2.63,  _ 

DATA  DVT2/2.  C7  ,2.U4,2-65  ,2 

DATA  DVT3/2.  07  ,2.Ut  ,  ^.67,  " 

DATA  " 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 
C 
C 

C CALL  DISSELA  ECUTINES  FOR  PLOT 

C 

CALL  TEK6  18 

CALL  MEDEUF 

CALL  RESET  (3HALL) 

CALL  HWSCAL  (  'SCREEN') 

CALL  PAGE  (1^.0,  S.5) 

CALL  GEACE  (CO) 

CALL  PHYSCR  (1.0-1.2) 

CiLL  AEEA2D  (  1  C.  0  ,b.  5) 

CALL  FBAME 

CALL  SWISSL 

CALL  BASALF  ('I/CSTE') 

CALL  MIXALF  (' STAND' ) 

CALL  INTAXS 

CALL    SHDCER     (  .  90  ,  1 ,.  0  15  ,  1  ] 

CALL    HEIGHT     (.  16) 


CALL    HEIGHT     (.  lb) 

CALLXNAME     ('(E)ISK     (LiOAEINGl  •  .  1  00] 

CALL    YNAKE     (  '  (  C  [  CEr  FICENT    OF     (1)  IFTS  '  ,  1  00  ) 

CALL    HEIGHT     (.290) 

CALL    MESSAG     ('  (H)  ELI  COPT E R     (C)ARPET     (P)LGTSS', 
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1  ICO. 2. 25, 7. 55) 

CALL      IlLXUGl         (./()) 

L.ui.     'j  r  d  r      l  _  .  L  ,  .   1  ,  -  .  r  ,  .  _  ,  .  U  5  , 

I*GRID=5 

CALL    THKCBV     (.C  15) 

CALL    PARA3 

CALL    LZGL1N 

CALL    CURVE     (DI -C6CC , 10, 0) 

(DL,Cc25,  10,  0) 

DI  ,C65C, 10,0! 

DL,Col5,  10,0 

DL,C7CC,  10,0 


t) 


CALL  CURVE 
CALL  CURVE 
CALL  CURVE 
CALL  CURVE 
CALL  DASE 
CALL  CURVE 
CALL  CURVE 
CALL  CURVE 
CALL    CURVE 


DVT1  ,C13,5,0) 
DV12,CLR,5, 0 
DVT3,CL3,5,0 
DVTm,CLR ,5,0) 

CALL    CURVE     (D V  15  .  CIR ,5 , 0 ) 

CALL    THKCBV     (.  C30) 

CALL    DCNEEL 

STCP 

ENE 
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c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c**  * 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


VARIABLES: 


CLE 

VT 

DL 

■ 
BB 


EESIGK    MEAN    LIFT     COEFIICENT 
TIP    VELOCITY 

SEIGHT     IS    CALCULATED     EIC^    POKES    EQUATION 
CSEFUL    LOAD    ELUS     EMPT^    WEIGHT 


ECWER     AVAIIAELZ    I  N*  HC  RS  EFOWE  E 


CL3  EQUALS  THE 
CL4  EQUALS  TbE 
DL5  EQUALS  THE 
DL6  EQUALS  TEE 
EL7  EQUALS  THE 
S'3  EQUALS  THE 
K<4  EQUALS  THE 
H5  EQUALS  THE 
B6  EQUALS  THE 
W7  EQUALS  THE 
WVI1  EQUALS 
WVI2  EQUALS 
WVI3  EQUALS 
WVTu  EQUALS 
•  VT5  EQUALS 
CVT1  EQUALS 
DVI2  EQUALS 
EVT3  EQUALS 
DVT4  EQUALS 
EQUALS 


DISK 
DISK 
DISK 
CISK 
DISK 


LOADING 
LOADING 
LOADING 
LOADING 
LOADING 


FCF 

FC5 
FCF 
F.CF 
FCF 


CLR=.3 

CLR=.4 
CLR=.5 
CLE=.o 
CLR=.7 


** 
** 
♦  * 
** 
***** 


SEIGHT    FOE    CLR  = 

SEIGHT    FOE     CLR=.4 

SEIGHT    FOE    CLR=.£ 

SEIGHT    FOR    CLR=-6 

SEIGHT    FOR     CLR=.7 
HEIGHTS     AT    V  7=6  0  0 
WEIGHTS    AT    VT=o25 
WEIGHTS    AT    VT=o50 
WEIGHTS    AT    VT=o75 
WEIGHTS     AT    VT=700 
TfcS    CCEHESPONCING    EISK    LOADING 
TEE    CORRESPONDING     EISK    LCALING 
TFZ    CORRESPONDING    CISK    LOADING 
TEE    CORRESPONDING    DISK    LOADING 
Er    CORRESPONDING    EISK    LCALING 


n  i. 

AT 
AT 
AT 


VT=600 
VT=o25 
VT=d50 
VT=o75 
VT=730 


EVT5    . 

REE    SOIOB    DISK    ECGNDAEY 

WBG    MAX     GFOSS     SEIGHT     BOUNDARY 
C************************** *********  ************************************ 

BEAL*4     DL3  (5)  x  EL4  i5)  ,DL5  (5)  , DL6J5)  fDL7  (5)  ,  »3(5)  .  S4  (5)  ,  ED3  (2)  , 
*«5  (5)  ,  Wo  (5)  ,  B  /  (5)  ,SVT1  (5)   ,WV~2  |5]  ,WVT3  (5)   ,  KVT4  (5)   .WtlG(z)  , 
*WVI5  (5)  ,  EVT1  (5)  ,DVT2  (5)  ,  D VI 3  (  E)  , DVT4  (  5)  ,D  VT5  (5  >  ,  L>  1  (  2)  ,  D2  (  2  ) 
C 
C DEFINE    DATA       

DATA    DL3/2. 06, 2. 07,2.07  ,,2.07, 2. 05/ 
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DATA  DI/2.52,2.9/ 

IJIU  Ntl£/2450.  ,2h:u./ 

UA'iA  i)Z/ 1.36  ,2.Z^/ 

DATA  RD3/2300.  ,2450./ 


C CALL    DISSELA     HCLTI1.ES    FOE    PLOT 

c 


1 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 

CALL 
CALL 
IIG2 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


618 
EUF 
ET 
CAL 
( 

CE 
SCE 
A2D 

me 

SL 
ALF 
ALF 
AXS 
CKS 
CKS 
CBH 
GET 
ME 
ME 
GHT 
DIN 

GRT 
i     ( 

CEV 

A3 

LIN' 

VE 

VE 

VE 

VE 

VE 

E 

VE 

Vr 

VE 

VE 

VE 

CEV 

DCT 

VE 

VE 

EEL 


(3HALL) 

(  'SCREEN') 
1 2 . C ,    9.5) 

1. 0,1.2) 


(0-0 


•  L/CSTC') 
'STA&D' ) 


(.90,  1,.015,  1 


16) 


('  (D)  ISK     (L)  OACIMGi  ',100) 

('   (G)SCSS     (W)  EIGHT     |  (LBS)  )  S' ,100) 

(.290) 

('  (H)  EIICOPTSB     (C)^EPET      (E)LOTS3', 
1CC,  1.0,1) 

(.20) 
2.  0,.  1,3.  C, 23  0 C.  ,5 C, 2550.) 

(-CIS) 


<DI3,K3,5,0) 
(DI«,H4,5,0) 
pL5,H5,5,0 

(DLc,«c,3  ,0) 
(  D  I  7  ,«/,  5  , 0) 

(DVT1  ,kVT  1,5,0) 
(DVT2,*VI2,5, 0) 
(DVT3,i'vT3,5,0) 
(DVT4,SVT«4,5,  0) 

(DVT5,'.VT5,5,  0) 
N.C35) 

(D1  ,WKG,2,0) 
(D2  ,3LE,2,0) 
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J;*********************************************************************** 

^♦♦♦••••..^^  GHAPHICAL     fcELICCPTcr    uuiui.  tfUbtiail                                 **„,*,,..,^ 

C*****»**                                   EAtlLY    OF     SOLUTIONS  ********* 

C********                                  CAEEET    PLOT    NUHEEH    5  ******** 

C********                                             21    AL    HANSEN  ******** 

r**»***»*  ******** 

..,.,,..  THIS    pbcGBAM    IS    DESIGNEE    TO    ILLUSTRATE    THE    FAMILY       ******** 

c********  c?    SOLUTIONS    FOE    HOVEF    AND     CSEFUL    LOAD                                ******** 

..,,..,.  BEQUIBIHIKSS    OF     A    TEETERING     FOTCR    SYSTEM                          ******** 

,,»,,.»,  HITH     pcxc?     DIAMETER,     FAX    GROSS     WEIGHT    AND                         ******** 

..,,.».,  ASEECT    RATIO     ECL'NDARIES.                                                         ******** 

C**»»**»*  CESEEVATION    CLASS    HELICCETEE                                           ******** 

C**********  ********** 

C******»*************»***** ***************  ****************************** 

C  **************************  **  *******♦*****.****************************** 

c**  **  ***** 

C**  NOMENCLATURE  *** 

C*  ** 

C*  VARIA3LES:  ** 

C*  ** 

C*  CLE  DESIGN    BEAK    LIFT     COEFIICENT  ** 

C*  VI  TIF    VELOCITY  ** 

C*  DL  DISK     LOADING 

C*  U  WEIGHT     AS    CALCULATED    FFCM    POKER    EQUATION  ** 

C*  HE  DSEFOL    LOAD    PLUS     EMPTi    HEIGHT  ** 

C*  PA  ECWER     AVAILABLE    IN    HCRSEtOWEE  ** 

C*  ** 

c****  ***** 

C*  EL3    EQUALS    THE    DISK     LOADING    FCF    CLR=.3 

C*  DL4    EQUALS    THE    DISK     LOADING    FCf    CLE=.4 

C*  DL5    EQUALS    TEE    DISK    LOADING    FCF    CLR=.5 

C*  DL6    EQUALS    THE    DISK    LOADING    FCF    CLR=.6 

C*  CL7    EQUALS    THE    DISK    LOADING    FCF    CLR=.7 

C*  H3    EQUALS    THE     WEIGHT    FOE    CLR=.3 

C*  K4    EQUALS    THE     "WEIGHT    FOR    CLR=.4 

C*  H5    EQUALS    THE     WEIGHT    FOR    CLR=.5 

C*  "W6    EQUALS    THE     WEIGHT    FOE    CLR=.c 

C*  H7    EQUALS    THE     WEIGHT    FOR    CLR=-7 

C*  HVT1    EQUALS    HEIGHTS    AT    VT=600 

C*  WVT2    EQUALS    WEIGHTS    AT    VT=q25 

C*  WVT3    EQUALS    WEIGHTS    AT    VT=650 

C*  WVI4    EQUALS    "WEIGHTS     AT    VT=675 

C»  HVT5    EQUALS    WEIGHTS    AT    VT=700 

C*  DVT1    EQUALS    IHE    CORRESPONDING    DISK    LOADING    AT    VT=600 

C*  DVT2    EQUALS    THE    CCRR  ESPO H  DI NG    DISK    LOADING    AT    VT=625 

C*  DVT3    EQUALS    THE    CORRESPONDING    DISK    LOADING    AT    V?=o50 

C*  EVT4    EQUALS    THE    CCRR ES? ON DI NG    DISK    LOADING    AT    VT=o75 

C*  DVT5    EQUALS    THE    CORRESPONDING    DISK    LOADING    AT    VT=7G0 

C*  EDE    EOTOE    DISK    ECUNEARY 

C*  WMG    MAX    GE0S5     WEIGHT    BOUNDARY 

C*  AR    HEIGHTS    FOE    THE    ASPECT    RATIC    30UNDARY 

C*  D3    DISK    LOADING    CCC  E  ES?  ON  DI  NG    TO    THE     ASPECT    RATIC    "WEIGHTS 

£***«***»***********»****»»********************************************* 

REAL*4    DL3(5),DL4(5),DL5(5)  ,DLt(5),DL7(5),W3(5).W4(5)  ,RDB (2)  , 
*i5  (5)  .  H6  15)  ,  W7  (5)  .WVT1  (5)  -HVT2  (5)  .WVT3  (5)  ,  KVT4  (d)     "WKG  (2)  , 
*HVT5  (5)  ,  EVT1  (5)  ,DVT2  (5)  ,  D  VT3  <  5)  ,  DVT4  (5)  ,  D  VT5  (5)  ,  D  1  (  2  )  ,  D2  (2  )  , 
*D3  <5)  ,  AR  (5) 
C 

C DEFINE    DATA       

DATA    DL3/2.06,2.07,2.07,2.07,2.05/ 

DATA    W 3/2  505-  ,2476. 2, 24  :>5. 46, 2  43  3. 85 ,24  18.25/ 

DATA    DL4/ 2. 42. 2.4^,2. 46, 2. 47, 2. 4 7/ 

DATA    W 4/2 467.  17,244  2.95,2  41 9.62,2399.72,2382-9  9/ 

DATA    DL5/2. 63, 2. 6 5, 2. 67, 2. o 9, 2.o9/ 

DATA    W5/2440. 72,24 1 9.51 ,2399. CC.237 9. 13,2 3o4. 65/ 

DATA    DL6/ 2.  74  .2. 7"i  .2.  79,  2. 81,  2. 81/ 

DATA    H6/24 24. 66, 2402.46,2 33 3. 56. 2 365. 10, 23 52./ 

DATA    DL7/ 2. 8 0,2. 6 3,2. 65, 2.8b,  2-87/ 

DATA    H 7/2411. 92,2390.72, 2 37 2. EC, 2 357-35, 23 42. 2 5/ 

DATA    HVT  1/25 0  5.  ,246  7.  17,  2  44  0. 7  2,2424.67,24  11.92/ 

DAIA    HVT 2/2476. 2,  24 4 2- 95.  24  19. 51 .2402.53,2390.72/ 

DATA    HVT3/24  55.4o,i4  19.  b2,239S.,2383.58,2372.S/ 

DATA    HVT 4/ 24 3 3. 8  5, 2 3 99. 7 2  ,2379.13,2365.  1,2357. 35/ 

DATA    W'VT  5/ 2416. 25, 2362.  99, 2364. 65, 2352. ,2342.2  5/ 

DATA    DVT  1/2.  06, 2. 42, 2. 63,  2.  74,2.6  0/ 


.,2U50./ 
3/ 


■ —    CALL    DISEELA     FCCTINES    FJR    PL01 


1 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
I 

CALL 
CALL 
IIGB 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STCP 
INC 


TEK618 

MECEGF 

RESET     (3HALL) 

HHSCAL     ('SCREEN') 

PAGE     (12. C,     S.5) 

-RACE     (O.C) 

(1.6,1.2) 


PHXSCR 
A  RZA2D 
FaACE 


(1  C. ,6.5) 


S»I££L 

BASALF     ('L/CSTE') 

.MIXALF     (•STAhC) 

INTAXS 


INTAXS 

Y1ICKS 

XIICKS 

SHDCHR 

HEIGET 

XNAKE 

YNAEE 

HEIGHT 

HEAI1N 


I?) 

j. 90,1,. 015,1) 

(.  1b) 
'   (D)  I 


SK      (L)  OAEING3',  100) 

CSS     (rf)  EIGHT     I  (LBS)  )  $«  ,100) 

) 


IOPTE5     (C)  iErET      (?)  LO: 


|    U    I     H    V.  0  .1  (Tl 

(.290) 

('  JH)  ELICOP 
ICC, 1.0,1) 

HEIGET     (.20) 

GRAF     (2.  C,.  1,3.  C,230C.  ,5C.,2550 

D  =  5 


'55 


GRAF 
ID  =  5 
THKCRV 
PARA3 
LEGLIN 
CURVE 
CURVE 
CURVE 
CURVE 
CURVE 
OASE 
CURVE 
CURVE 
CURVE 


(-C13) 

(DL3,V.3,5,0) 
[014,-4,5,0) 
DL5,«5,5,0) 
(DLb,'*c,5,0) 
[DL7,W~,5  ,0) 

(D\11,SVT  1,5,0] 

*»  V  T  2 ,  5  ,  0  j 

ti  rr    1        C  i 


CURVE     (DV12,'«VT2,5,0 
CURVE     (D \ 13, « VI 3,5,0) 
CURVE     (DVl4,'.VT4,5,  0) 
CURVE     (D\15.UVT5,5 , 0) 
THKCF.V     (.030) 
CHNCCT 

CURVE     ( C  1  ,  W  B  G  ,  2  ,  0) 
CURVE     (D2, REE, 2,0) 
CURVE    (D3,£H,5,0) 
DONEEL 
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C       TEIS    PROGRAM    IS    CESIGN    TO    GENERATE    THE    LATA    FOR    THE    GRAPHICS 

L       sttUUON     or     in     ariGn-i     ANn    thf     USEFUL    LOAD    EOUATIOH.        Taj"    tc     n.uc 

l_       fiflil    SXJSf    IN    A    LAbtiT    ELOT     HELICCETER    DESIGN    PARAMETRIC    OPTIMIZATION 

C 

C  ASSUMPTIONS: 

C  1>     ENGINESPECIFIEI 

C    VAHIAELE    OPTIONS 

C 

BEAL*4    ClE,PA,tI,K1,K2,K3,K4,i<5,E,S,A,E,;i<10),WB(10) 

INTEGER     VT, D,I,CI 
C 
c 

CALL    FRTCKS     ('FliEEEF    ','02  «,'EISK  ,,»CaPT1', 

>  'EATA  ',  •  A  ') 

C 

C 

CALL    FBTCES     ('FILEEEF    • , • 03                  ".'DISK             ,,'CEPT2,# 
>'EATA              ','A                       ») 
C 

c 

C  CLB=    DESIGN    MEAN    LIFT    COEFFICIENT 

DC    90    CL=3.7 
CLB=CL*  (C.  1) 
BRITE(2, 10) CLE 
C  EA=    POWEB    AVAILAELE     EP 

PA-206 
C  £1=    DISK    LOADING 

C  VT=    TIF    VELOCITY    FT/SEC 

C 

C  CONSTANTS    BASIC    ON    CIB 

C 

K1=  (0.  9563)  /CIS*  (1+1. 8078  *CLB**2) 
K2=PA*10**5/K 1 

K3=(G.00C25929)/CLB»  (1+  1.8078*CL3**2) 
KU=55348C.0/K1 
K5=3695.7/K1 
DC    100       D=200,300 
EL=D*  (0.C1) 
1  =  0 

DO    11C      VT  =  6CC  ,70C,25 
C 
C  ARRAY    INCRECENTER 

1=1  +  1 
C  WEIGHT    EQUATION 

C 

W(I)  =  (K2*(1-(411.E1*EL**.75)/(\T**1.5)*(1  +  K3*VT/DL**„5)**.5)-K4) 
1/  (VT*K5*EI**.5) 
C 

C  CALCULATION    OF    WE    LATA 

C 

A=»(I)  /DL 
B=  (A/3.  14) **. 5 
P=A**. 5/V1 

S= (6.* EL)/ (0. 0 023679 »CLR  *VT**2) 
C 

C  ASSURING    A    TEEIESING    SYSTEM 

C 

HE  (I)  =17  17. 9+ 122  1.*?**.  8  63+26€.*P**. 7+1744  9./VT** 1.  14 
1  +  2.886*A**.5*P*».57+.  191  *'»(!)  **.916<-.0294*«(I)  **.  9  9 
2*35.  15*  (W  |I)/VT)  *  JA**.1  85) *S**. 33  +  .  0Q88*S  (I)  *A**.21 
C  W8ITE     (5, 2C) VT 

C  WRITE  lb, 30)  S 

110  CONTINUE 

WRITE     (2,3  1)  Dl, '»  (1),  WE  (1  )  ,t>  (2)  ,W  E  (2) 
WRITE     (3,:1)     EI,  k  (3)  ,  «3(3)  ,'i  (4)  ,  »  E  (4)  ,  W  (5)  ,WB  (5) 
100       CONTINUE 
90       CONTINUE 
C  FCBMAT    STATEMENTS 

10  FOEMAT     ( ■  1'  , ' CLB=' ,  1F10.4///) 

31  FORMAT     (6  (2X,F  10.3)  ) 

STOP 
ENE 
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C         ZE2*i      cavavBD     .lo      Lbit.i      xvj     uuIAxh      iaL     tiitii^     <:iun.      r  H  -J  0     TrlL  * 

C       ASPECT    BATIC    GSAPfc'    INTERSECTION     ECINT    DATA.    IN    ORDER    TO    SBAPH  * 

C       I  HE    ASPECT    HATIC    ECUNEABY    ON    THE    CAIN    CARPET    PLOT    A 


C       CCBBESFCNCINC    WEIGHT    BEST    BE    OBTilNEB.     THIS    FBOGHA2!    UTILIZES  * 

C       THE    iiEIGHT    ESTIMATION    FCHMOLATIOh    FBOX    THE    CARPET    PHOGRAM.  * 

C************»*»»***^»**»**»*»*****»******  »«**#*»*****»<■*«»********»* 
C ****#»**#*****»***#***»»♦* ****»***^*****#*»**»*****************#**** 
C       ASSUMPTIONS:  * 

C  1>    ENGIRESPECIFIEI  * 

C    VABIAELE    OPTICNS 

BEAL*4     VI  (5)  -CIJ5)  ,CIB{5)  ,W  (5) 
DATA    CLR/. 625,. 57 5,. 5 35,  .46  0,.  4 2/ 
DATA    VT/6C0.  -625-  ,6  5 C.  ,67  5.  , 7CC./ 
DATA    DL/ 2. 75, 2. 7 5, 2. 72, 2. 67, 2. 53/ 

CALL    FBTCHS     ('EILEBEE    ','03                  ','DISK             ^'CEPTAS1, 
>«DATA             •  ,'A                     ') 
C 

c 

DC    90    L=1.5 

SRITE{3,  1C1CLP  (L) 

HEITE(3,2C    VT  (I) 

¥BITE(3,3C)  DL  jlj 
C  PA=    POKES    AVAILABLE    HP 

PA=206 
C  CL  =    DISK    LOADING 

C  Vl=    TIP     VELOCITY    ET/SEC 

C 

C  CONSTANTS    BASEL    CN    CIK 

C 
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31 
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END 

STBTEHEKTS 

('1'     •CLR=,,1F10.U///) 
2X,f VT  =  «  ,  1F1C.4///) 
2i,«DI='  ,  1F1C.4/// 

(2X,F1C3) 
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APPENDIX  D.   PROGRAMS  TO  ACCESS  HESCOMP 

This  section  contains  the  control  language  programs 
needed  to  access  HESCOMP  on  the  IBM  main- frame  computer 
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APPENDIX  E:   HESCOMP  INPUT  AND  OUTPUT  EXAMPLES 

This  section  contains  samples  of  the  IBM  computer 
input  and  output. 
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